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We show how dipolar interactions between dysprosium atoms in an optical lattice can be used to 
obtain fractional quantum Hall states. In our approach, dysprosium atoms are trapped one atom 
per site in a deep optical lattice with negligible tunneling. Microwave and spatially dependent 
optical dressing fields are used to define an effective spin-1/2 or spin-1 degree of freedom in each 
atom. Thinking of spin-1/2 particles as hardcore bosons, dipole-dipole interactions give rise to boson 
hopping, topological flat bands with Chern number 1, and the v — 1/2 Laughlin state. Thinking 
of spin-1 particles as two-component hardcore bosons, dipole-dipole interactions again give rise to 
boson hopping, topological flat bands with Chern number 2, and the bilayer Halperin (2,2,1) state. 
By adjusting the optical fields, we find a phase diagram, in which the (2,2,1) state competes with 
superfluidity. Generalizations to solid-state magnetic dipoles are discussed. 

PACS numbers: 67.85.-d, 73.43.-f, 37.10.Jk, 05.30.Pr 


I. INTRODUCTION 


In addition to their fundamental importance, topolog¬ 
ical phases of matter may eventually enable the realiza¬ 
tion of fault-tolerant quantum computing [1] and robust 
quantum state transfer [2]. Thanks to unprecedented 
controllability and purity, synthetic atomic, molecular, 
and optical systems are rapidly gaining momentum in 
their ability to exhibit some of the most exotic topolog¬ 
ical phases [3-10]. Dipolar systems, such as Rydberg 
atoms, polar molecules, and magnetic atoms, have re¬ 
cently attracted a particular degree of attention thanks 
to the strength of the interactions and the natural link 
between dipolar interactions and topology [11-18]. In¬ 
deed, thanks to the Einstein-de-Haas effect, dipolar in¬ 
teraction can convert internal angular momentum into 
angular momentum describing the rotation of the two 
interacting dipoles around each other. It is not surpris¬ 
ing that the resulting rotation can be harnessed for gen¬ 
erating effective gauge fields. However, nearly all dipo¬ 
lar topological literature to date focuses on interacting 
topological states arising from flat Chern bands with a 
Chern number equal to one. The exception is Ref. [17], 
which does consider bands with higher Chern number, 
but which provides no path for making the bands suffi¬ 
ciently flat for realizing interacting topological phases. In 
this article, utilizing the magnetic atom dysprosium [19], 
we demonstrate for the first time how to create a topo¬ 
logical flat band with Chern number 2 and to utilize the 
resulting band-structure to realize the bilayer Halperin 
(2,2,1) fractional quantum Hall state [20]. By adjusting 
the magnetic field, we find a phase diagram, in which 


the (2,2,1) state competes with neighboring superfluids. 
More generally, our goal here is to demonstrate the con¬ 
trollability and potential of magnetic dipoles for simulat¬ 
ing many-body phases; however, we emphasize that the 
constraints required by our approach are challenging for 
current generation experiments, and rather than an ex¬ 
perimental blueprint, our approach should be viewed as 
a general framework. 

The remainder of the article is organized as follows. 
In section II, we describe the effective Hamiltonian as¬ 
sociated with a two-dimensional lattice of ultracold dys¬ 
prosium atoms. We demonstrate that magnetic dipo¬ 
lar interactions mediate both long-range dynamics (hop¬ 
ping) and interactions. Microwave and optical radiation 
is used to break time-reversal symmetry and to con¬ 
trol the specific nature of the atomic degrees of free¬ 
dom. By tuning these dressing parameters, we realize 
the v = 1/2 fractional Chern insulator and clarify its 
characteristics with a variety of numerical diagnostics. 
In section III, we generalize our approach to solid-state 
magnetic dipoles. Specifically, we consider the example of 
Nitrogen-Vacancy defects in diamond and offer a route to 
sub-optical-wavelength resolution dressing via patterned 
dielectrics. In section IV, we discuss how to realize a 
C = 2 topological flat band with dysprosium, by consid¬ 
ering an effective spin-1 atomic degree of freedom. Upon 
populating this band-structure with a finite density of in¬ 
teracting particles, we find a ground state that exhibits 
a Hall conductivity a xy = 2/3 consistent with the (2,2,1) 
Halperin state. In section V, we elaborate on the exper¬ 
imental considerations and discuss the challenges using 
current technologies. Finally, in section VI, we present a 
brief outlook. 
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FIG. 1. (a) Dysprosium atoms are loaded one atom per site 
into a square lattice in the X-Y plane, (b) The quantization 
axis of the atoms z (red) is determined by the polarizations 
of applied fields. To obtain the x-y-z coordinate system, one 
rotates the X-Y-Z coordinate system by <E>o around the Z 
axis and then by ©o around the y axis, (c) The relevant 
161 Dy level structure used for obtaining the flat band with 
Chern number C — 1 and the corresponding v — 1/2 Laugh- 
lin state. The lower hyperfine levels (nuclear spin I = 5/2) 
are in the J — 8 ground electronic state. Far-detuned light 
(not shown) provides a tensor AC Stark shift oc M 2 with an 
F-dependent coefficient. The 421 nm light (purple), which 
drives a transition to the F' s= 17/2 hyperfine level of the 
4/ 10 ( 5 /8)6s6p( 1 Pi°)(8, l)g state, and the microwaves (black) 
define the dark dressed states |f) and |f). The dressing opti¬ 
cal field is different on the two sublattices [green and yellow 
in (a)] making the state |t) sublattice-dependent. 


II. C = 1 FLAT BAND AND THE v = 1/2 
LAUGHLIN STATE 


To introduce the main features of dysprosium, as well 
as dipolar-mediated topological flat bands and fractional 
Chern insulators, we will first show how to obtain a flat 
band with Chern number (7 = 1 and use it to realize a 
v mm 1/2 Laughlin state with dysprosium. As shown in 
Fig. 1(a), we consider loading 161 Dy one atom per site in 
a square lattice in the X-Y plane with nearest-neighbor 
spacing of Ai at /2 = 266 nm created with off-resonant 
light of Ai a t = 532 nm wavelength. We further assume 
that the lattice is so deep that tunneling is negligible. As 
shown in Fig. 1(b), the quantization axis £, determined 
by the polarizations of applied fields, is pointing in a tun¬ 
able direction different from Z. An off-resonant linearly 
polarized optical field (not shown) provides a tensor shift 
oc M 2 shown in Fig. 1(c). Optical and microwave dress¬ 
ing fields are then used to construct dressed dark states 
1 1) and | 7 ), where the site-dependent coefficients 8 , v, 
and w in the definition of |t) are controlled by the am¬ 
plitudes of the applied fields. We choose to work with 
M near zero to maximize the transition dipole moments 
between the states involved and to simultaneously mini¬ 


mize the sensitivity of the energy levels to stray magnetic 
fields. 

Interactions between the effective spin- 1/2 particles 
are mediated by magnetic dipole-dipole interactions. 
Ignoring the small nuclear dipole moment, the inter¬ 
action between two dysprosium atoms i and j sepa¬ 
rated by (R, #, </>) in the spherical coordinates associ¬ 
ated with the x-y-z coordinate system, in the units of 
Fo{9jFb) 2 /( 47 tR 3 ) (where the Lande g-factor is gj = 
1.24 [21]), is given by 


Hij = (1 — 3 cos 2 6) 



jFj - + FF\ 

2i + J+ J+ + h.c. . 


( 1 ) 


Here J = 8 is the total electronic angular momentum 
associated with the electronic ground state. We set h = 1 
throughout the paper. Assuming that the dark states 
|t) and |t) are separated from the bright states by an 
energy D^ r /A (here Odr and A being the Rabi frequency 
an detuning of the optical dressing fields; the microwave 
Rabi frequency is assumed to be ~ H^ r /A) that is larger 
than the dipole-dipole interaction strength, we project 
the Hamiltonian onto the four-dimensional Hilbert space, 
where each of the two atoms is in state |t) = | 0 ) or |t) = 
a |0). We have introduced the creation operator at for 
a hardcore boson on site j. The Hamiltonian then takes 
the form 


H ( 


dd 




( 2 ) 


where hi = dja* and where the hopping amplitude Uj 
and the density-density interaction Vij (see Appendix 
A) depend on the vector (R, #, (j )) separating atoms i 
and j, the direction ( 0 q,$o) °f the quantization axis 
[see Fig. 1(b)], the dressing coefficients Sij, Wi,j 
[see Fig. 1(c)], permanent dipole moment /14 = (4|J 2 |4) 
and transition dipole moments ^ 53 , // 13 , and [126 (here 
hag — (ct\J z,± |/3), where the difference between the M 
quantum numbers of |a) and \/3) determines the choice 
of z or ±). Notice that, in contrast to electric-dipole im¬ 
plementations, the bare | F,M) states in Fig. 1(c) have 
substantial dipole moments even in the absence of ap¬ 
plied fields. 

In Fig. 2, we show the (7 = 1 topological flat band 
and the resulting v = 1/2 Laughlin state on a torus, ex¬ 
hibiting the expected gapped two-fold degenerate ground 
state. Each ground state was verified to have the many- 
body Chern number of 1/2 and the correct quasi-hole 
statistics obeying the generalized Pauli principle [ 22 ]. 
The specific values, used in Fig. 2, of site-dependent co¬ 
efficients 8 , v, w and of the direction of (@o,ffio) of the 
quantization axis relative to the X-Y-Z plane are given 
in Appendix A. In particular, we would like to set arbi¬ 
trary s, v, and w on one sublattice [yellow in Fig. 1(a)] 
and arbitrary 8 , v, and w on the other sublattice [green 
in Fig. 1(a)]. On top of that, we would like the value of 
w to alternate every other row. The 421 nm optical fields 
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FIG. 2. (a) Flat topological band with Chern number C — 1. 
The flatness of the band (band gap divided by band width) 
is « 8. (b) The momentum-resolved eigenvalues for the 

v = 1/2 fractional Chern insulator with 6 hardcore bosons 
on a 6 X 4 torus. The momentum sector n 2 + 4ni corre¬ 
sponds to ( k x ,k y ) = (ni/3, 712/2 — ni/3)7r, n 1 = 0,1? 2, and 
ri 2 = 0,1,2,3. The spectrum features a gap separating two 
degenerate ground states at (k x ,k y ) = (0,0) and (0,7r) (blue) 
from the other states by a gap. 


shown in Fig. 1(c) are key to generating this site depen¬ 
dence of 8, v , and re, while the microwave fields provide 
spatially uniform couplings. Specifically, first, we apply 
an optical field with a nonzero 7r-polarized component 
(with respect to the quantization axis z) that is uniform 
across the square lattice by having the light’s ^-vector 
perpendicular to the plane of the atoms. Furthermore, 
the beam’s polarization can always be chosen such that 
its cr + and a~ components (with respect to the quan¬ 
tization axis z) have equal intensity. These <j ± compo¬ 
nents will thus only result in AC Stark shifts that keep 
the energies of |F, M) and | F,—M) equal and that can 
therefore be simply absorbed into the tensor shifts. Sec¬ 
ond, since 421 nm is less than the wavelength of light 
Aiat = 532 nm used to create the optical lattice, we can 
take a beam propagating along X and rotate its ^-vector 
slightly around Y, such that, within the XY plane of the 
atoms, it acquires periodicity of twice the lattice spacing 
along X. We will have two such beams corresponding to 
clockwise and counterclockwise rotations around Y. We 
then similarly take a beam propagating along Y and ro¬ 
tate its k- vector slightly around X, such that, within the 
XY plane of the atoms, it acquires periodicity of twice 
the lattice spacing along Y. Again, we will have two such 
beams corresponding to clockwise and counterclockwise 
rotations around X. The electric field of each of the re¬ 
sulting four beams is described by two complex numbers 
corresponding to the amplitudes of the two transverse 
polarization components. The resulting eight complex 
numbers can be generically tuned to get no 7r-polarization 
amplitude on any of the atoms and simultaneously arbi¬ 
trary cr + and a - amplitudes on the two sublattices sub¬ 
ject to the desired sign alteration every other row. 
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FIG. 3. (a) Magnetic field required to tune the hyperhne 

coupled NV states to their desired resonances, (b) Optical 
dressing M-scheme which enables sufficient control to realize 
topological flat bands. 

III. IMPLEMENTATION IN NV CENTERS 

To depict the generality of our construction, we 
now consider magnetic dipoles associated with Nitrogen- 
Vacancy color centers in diamond. The NV center has 
received a tremendous amount of interest in recent years 
owing to the fact that its electronic spin can be polarized, 
manipulated and optically detected under ambient con¬ 
ditions [23-29]. Each NV center also harbors a localized 
nuclear spin, which exhibits extremely long coherence 
times [30]. Using a combination of these electronic and 
nuclear degrees of freedom, we will demonstrate the abil¬ 
ity to realize topological flat bands. Our approach will 
be analogous to the previous section: namely, the use of 
microwave and optical fields to break time-reversal sym¬ 
metry and to realize appropriately dressed eigenstates. 

The electronic ground state of each NV center is a 
spin-1 triplet described by the Hamiltonian, 

H nv = D 0 Sz + ti e BS z , (3) 

where Dq/2tt = 2.87 GHz is the zero field splitting, fi e = 
— (27r)2.8 MHz/G is the electron spin gyromagnetic ratio, 
and B is a magnetic field applied parallel to the NV axis. 
This electronic spin is coupled via hyperfine interactions 
to the 1 = 1/2 nuclear spin of the 15 N impurity via 

Hhf = A\\Sj z + A ± (S x i x + Syly), (4) 

where A\\/2ir ~ 3.0 MHz and Aj_/2tt ~ 3.7 MHz [31- 
33]. We assume that a dc magnetic field tunes the en¬ 
ergies of states |0, — |) and |l,^) to be nearly equal 
and simultaneously far-detunes the energies of states 
| — 1, db — ) , where states are labeled by | S Z ,I Z ). The A± 
term in Eq. (4) mixes the |0, |) and 11, — |) states, 
yielding the energy levels shown versus magnetic field 
in Fig. 3(a), where we have defined the eigenstates | A) = 

- « |°, I), | B) = | 0 , —1>, I C) = |l, §), and 

| D) =55 a 11, — + P 10, |). To allow for resonant hops 
of spin excitations, we work at B ^ 1028 G where states 
| B) and | C) are nearly degenerate. 

The effective states we use on each NV center are |0) = 

| A) and |1) = s\B) +v\C) +w\D). Analogous to the case 
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FIG. 4. Control of solid-state spins, (a) Two dimensional ar¬ 
ray of implanted NV centers near diamond 110 surface, with 
NV lattice spacing a — 20 nm. A uniform applied electric field 
and a dielectric nanostructure patterned on the surface modu¬ 
late the displacement field at each NV center, (b) Schematic 
optical dressing scheme. The left panel corresponds to the 
bare level structure and the right panel corresponds to the 
shifted level structure in the presence of gating. Spin states 
|q=l) in the electronic ground state are coupled to the | A 2 ) 
excited state with cf ± polarized light. We have left out the hy- 
perfine structure for simplicity. Modulated displacement field 
brings the two frequencies (red and blue) of optical dressing 
in and out of resonance, (c) Patterned dielectric on the sur¬ 
face of diamond (with an indium tin oxide back gate) which 
can yield a local dc field on the order of 10 3 V/cm as shown 
in (d). 


of dysprosium, the coefficients s,v,w are determined via 
an optical “M” dressing scheme [Fig. 3(b)] where the two 
excited states are |±) oc \E X )±\A2), with \E X ) , IA2) being 
two specific electronic excited states of the NV [34, 35]. 
The state |1) is the so-called dark state of the M-scheme 
with s = ^ 2 ^ 4 /^, v = w = —Q1Q4/Q. Note 

that lasers 1 and 3 must be linearly polarized, while lasers 
2 and 4 are circularly polarized. This elliptical polariza¬ 
tion of light explicitly breaks time-reversal symmetry. 

The mixing angle tan(^) = \si/vi\ characterizes the 
strength of the effective dipole moment of 11), thereby de¬ 
termining the magnitude of the interactions. In the limit 
0i 0, s 0; since the state \B) carries no electronic 
spin dipole moment, the dipolar interaction strength in¬ 
creases as 0i —>• 0. Topological flat-bands are found for a 
variety of parameter regimes, and fractional Chern insu¬ 
lating ground states are typically obtained for 6^ > 0.5, 
where long-range interactions are relatively weak. 


A. Optical Dressing and Strain 

One challenge that arises in the context of implement¬ 
ing topological phases with solid-state magnetic dipoles is 
the ability to vary dressed states on length scales smaller 
than an optical wavelength. Indeed, to generate to ob¬ 
tain flat topological band structures, we require spatially 
inhomogeneous optical dressing from site to site, a non¬ 
trivial task for lattice spacings a ~ 20 nm (required to 


achieve sufficiently strong magnetic dipole-dipole interac¬ 
tions) well below optical resolution. In the case of NVs, 
one can, in principle, accomplish this task by modulat¬ 
ing an applied dc electric field from site to site, using 
a patterned conducting nanostructure on the surface of 
the diamond as shown in Fig. 4. A dc electric field ap¬ 
plied parallel to the NV axis shifts the electronic excited 
state with respect to the ground state triplet due to its 
strong electric dipole moment, d^ s ~ (27r)l MHz cm/V. 
The ground state dipole moment is significantly weaker, 
d^ s ~ (27r)0.35Hz cm/V, and can be safely neglected 
[35]. Thanks to the application of a dc voltage between 
the patterned surface conductor and a back gate, the 
electric field on red and blue sites tunes the optical tran¬ 
sitions of the NV centers in and out of resonance with 
red and blue optical driving lasers, which are applied 
globally [Fig. 4(b)]. For a diamond sample of thickness 
10 fi m, an applied voltage of V ~ 1 V generates a dc 
field of E ~ 10 3 V/cm [Fig. 4(d), field difference between 
two x’s], resulting in excited state shifts of order ~ GHz, 
which is significantly larger than the intrinsic line width 
7e/27T ~ 10 MHz. 

While this approach, in principle, allows for a site- 
dependent optical dressing, one needs to be extremely 
careful to consider the effects of local strain fields. These 
strain fields couples to NVs in the same way as an ef¬ 
fective electric field. In diamond grown by chemical va¬ 
por deposition (CVD) and implanted with NV centers, 
unwanted defects unavoidably lead to local variations in 
strain. The two predominant sources of strain in CVD di¬ 
amond are point defects, such as vacancy clusters and in¬ 
terstitials, and line defects, such as stacking faults which 
align along the growth direction. 

To estimate the strain from point defects, we take a 
defect concentration of 1 ppm [36-38], corresponding to 
an average point defect separation of r pt ~ 20 nm. For 
a defect strength of A ~ 10 -4 nm 3 , the resulting strain 
variations are of order e^A/r 3 t ^ 10 —8 , inducing shifts 
of the NV excited state of ~ 10 MHz [39-41]. 

Next we turn to the expected dominant source of lo¬ 
cal strain variation, arising from stacking faults aligned 
along the diamond growth direction. We model the strain 
from stacking faults using elasticity theory and assuming 
straight line defects. For a typical Burgers vector of mag¬ 
nitude b ~ 2 A, the strain at distance fii n from the line 
defect is given by e ^ bfi/ [27r(l — z/)rn n ], where /i is the 
shear modulus in diamond and v is Poisson’s ratio [42]. 
For an areal density of 10 4 /cm 2 [43], corresponding to an 
average distance of ^ 100 /im between line defects. We 
estimate typical strain variations of 1-10 GHz on length 
scales of 100 /im, in agreement with recent experiments 
[44]. However, we expect much smaller strain variations 
on length scales smaller than r\in- In particular, in a 
250 x 250 nm region between line defects, numerical sim¬ 
ulations give average strain variations of ~ 10 MHz. 

These estimates suggest that local strain variations in 
a small NV lattice might not destroy the spatially inho¬ 
mogeneous optical dressing obtained via patterned elec- 
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trie fields; however, such an implementation is extremely 
challenging and our discussion of NV centers is meant 
mainly as a proof-of-principle analysis. 


IV. C = 2 FLAT BAND AND THE HALPERIN 
(2,2,1) STATE 
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Having demonstrated the generality of our approach by 
extending it to NV centers in diamond, we now switch 
back to the dysprosium implementation. Having intro¬ 
duced in Sec. II the potential of dysprosium for creat¬ 
ing fractional quantum Hall states on the example of the 
v — 1/2 Laughlin state and the underlying C — 1 flat 
band, we now move on to the construction of the topo¬ 
logical flat band with Chern number (7 = 2. At z/ = 1/3 
filling fraction, this band will give rise to the Halperin 
(2,2,1) bilayer fractional quantum Hall state. 

We will follow the idea of Ref. [45] for generating the 
C — 2 flat topological band. Suppose one has created a 
(7 = 1 flat topological band on a square lattice with a 
two-site unit cell for a boson a, like we did above for the 
hardcore bosons arising from the spin-1/2 model in dys¬ 
prosium. Now let us introduce another species b of bosons 
that obeys exactly the same hopping Hamiltonian, except 
it is shifted relative to the hopping Hamiltonian for a by 
one lattice site in the X direction. Diagonalizing the re¬ 
sulting Hamiltonian clearly gives four bands: two (7 = 1 
bands and two C = — 1 bands. The idea is then to merge 
the two flat (7 = 1 bands into a single flat (7 = 2 band. 
To do this, on one sublattice, one defines the hardcore 
boson A = a and the hardcore boson B = 6, while on the 
other sublattice one does the opposite and defines A = b 
and B = a. The resulting model recovers the full transla¬ 
tional symmetry of the lattice and has therefore only two 
bands with C = ±2. We will realize the vacuum state 
and the two species of bosons on each site using a spin-1 
particle. Therefore, our bosons will be hardcore both to 
themselves and to each other, interactions that will be 
sufficient for realizing the Halperin (2,2,1) bilayer state. 

We start with the same equation (1) as in the (7 = 1 
discussion, but now project onto dressed states ||), |j"), 
and |ff), whose precise construction in terms of the 
ground hyperfine states of 161 Dy is relegated to Appendix 
B. We then define |0) = ||), a 1 |0) = |t), and P |0) = |ff). 
Choosing the dressed states in such a way that long-range 
density-density interactions vanish and such that a) and 
ZF have the same hopping matrix elements t[- (corre¬ 
sponding to a C = 1 flat band) but shifted relative to 
each other by one unit in the X direction, we obtain 

HdA = (5) 

Following the above-described redefinition from (a, b) to 
(A, T>), we arrive at a Hamiltonian describing a C = 2 
flat band. In Fig. 5(a), we show the resulting (7 = 2 
flat topological band. At v = 1/3 filling fraction, we 


FIG. 5. (a) Flat topological band with Chern number (7 = 2. 
The flatness of the band (band gap divided by band width) 
is « 11. (b) The momentum-resolved eigenvalues for the 

v — 1/3 fractional Chern insulator with 4 two-component 
hardcore bosons on a 3 x 4 toms. The momentum sector 
corresponds to (k x ,k y ) = (2ni/3, 712 / 2 ) 77 , n± — 0,1,2, and 
n 2 = 0,1,2,3. The spectrum features a gap separating the 
three degenerate ground states at (k x ,k y ) = (0,0), (27t/ 3,0), 
and (47 t/ 3,0) (blue) from the other states by a gap. 


will show that this topological flat band gives rise to the 
bilayer Halperin (2,2,1) state [20] shown in Fig. 5(b). To 
understand the nature of this state, it is helpful turn to 
the usual A-matrix description of bilayer quantum Hall 


systems, with wavefunction, ]Ji<j( z u ~ z ijY Hi<j( z 2i ~ 

^2j) m Hi,j( z n ~ Z2j) n e~i^^ Zli ^ + ^' l* 2j 'l ^ where z pi is 
the complex lattice coordinate of the ith hardcore boson 
in the pth layer. Here, l and m are even integers so that 
the wavefunction is consistent with Bose statistics, while 
n can be any non-negative integer. The A-matrix is then 
defined as 



In our case, from the structure factor, one finds that the 
state does not seem to preference a specific orbital type, 
suggesting a charge vector q = (1,1). The transverse 
Hall conductance a xy is then given by 

cr xy = qK~ 1 q. (6) 

From the numerics, we observe a gapped three-fold de¬ 
generate ground state on a torus, each exhibiting a trans¬ 
verse Hall conductance, a xy = 2/3 [45, 46]. Combined 
with the fact that ground state degeneracy is given by 
the determinant of A, this suggests that we are indeed 
observing the (Z,m, n) = (2,2,1) state. 

In Fig. 6, we show the entire (7 = 2 phase diagram [as 
a function of electric field tilt angle (0o, 4>o)], where the 
state in Fig. 5(b) is a single point (shown as a red cir¬ 
cle) at ©o = 0.68 and 4>o = 5.83. Similarly to Fig. 5, the 
phase diagram is obtained using exact diagonalization on 
a 3 x 4 torus, with two orbitals per site, and a total of 
4 particles, corresponding to a density of v = 1/3 parti¬ 
cles per unit cell. The finite momentum superfluid has 
states which twist into the continuum but the ground 
state looks like it is at finite momentum instead of in 
the (k x ,k y ) = (0,0) sector. The phase diagram is (a 
schematic smoothed version) based upon a back grid of 
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FIG. 6. The phase diagram obtained by starting with the 
Halperin (2,2,1) state shown in Fig. 5 (red circle) and varying 
the direction (©o, <ho) of the quantization axis. 


100 points. We sample at 10 equally spaced points along 
©o and $ 0 and at each, we obtain the spectrum, the 
structure factor of the lowest energy state in all momen¬ 
tum sectors, the superfluid response from boundary con¬ 
dition twists, and the manybody Chern number a xy . In 
the region called (2,2,1), a xy is numerically 2/3 within 
< 1% error, which is a strong indicator that this is in¬ 
deed the Halperin (2,2,1) state. 

V. EXPERIMENTAL CONSIDERATIONS 

Our scheme relies on the following ladder of energy 
scales: = 1 GHz) >> (Ht/27r 300 Hz) 

((^ r /A)/27r = 30 Hz) > (y dd /27r = 3 Hz) > 
lE FQU /27r = 500 mHz) > (r sc /27T = r' c /27r = 
100 mHz). The first inequality Vhf V t ensures that 
the tensor AC Stark shift V t does not mix different F 
levels separated by Vhf, which, in turn, ensures that the 
dipole moments of the involved states and transitions are 
given by their zero-field values, which simplifies the cal¬ 
culations. The second inequality V t 0^ r /A ensures 
that the optical fields (Rabi frequency fJdr and effective 
two-photon Rabi frequency f^ r /A) and microwave fields 
(Rabi frequency f^ r /A) are sufficiently weak to spec¬ 
troscopically resolve the different transitions split by the 
tensor shift, allowing us to consider only the desired cou¬ 
plings. The inequality ^ dr /A Vdd ensures that dipole- 
dipole interactions Vdd, whose strength can be read out 
from the vertical energy scale in Figs. 2(a) and 5(a), do 
not perturb the dark states |t) and |^) (and |ff) for the 
C = 2 example) defined by the dressing fields and do 
not cause transitions from these dark states to the bright 
states. We will discuss below that the small 3 Hz in¬ 
teraction energy scale is not fundamental and can likely 
be significantly increased. The inequality Vdd ^fqh, 
where E’fqh is the fractional quantum Hall energy scale, 
which can be read out from the gap in Figs. 2(b) and 5(b), 
is not required but arises naturally since Efqh is deter¬ 


mined by Vdd, so we included this inequality in the ladder 
of energy scales. Finally, the inequality E’fqh r sc ,r' c 
ensures that the photon scattering rates r sc and T' sc due 
to the optical fields (used for and for Vt, respec¬ 
tively) are much smaller than the energy scale Efqh °f 
the Hamiltonian of interest. 

Let us now discuss the scattering rate r' c arising due 
to the optical fields creating Vt- In alkali atoms, the ra¬ 
tio of the tensor Stark shift Vt to the scattering rate is 
small, and, for the case of Cs, maximized at only 8.2 (a 
fraction of the ratio between the hyperfine splitting and 
the linewidth), when the light is tuned between the D\ 
hyperfine transitions [47]. On the other hand, the large 
orbital and spin angular momentum in the ground state 
of dysprosium gives rise to qualitatively stronger tensor 
light shifts and, hence, to much more favorable ratios of 
Vt to the scattering rate [48]. In particular, the tensor 
shift is on the order of Vt = If' 2 /A', while the decay 
rate is on the order of T' c = T f (Q f / A') 2 . Here 12', A', 
and T' are the corresponding Rabi frequency, detuning, 
and spontaneous emission rate. Assuming we use 741 nm 
light, r'/27T = 1.8 kHz [49]. To get r' c /27r = 100 mHz 
and Vt/27ir = 300 Hz, we need A' /2tt = 5 MHz and 
fy/27r = 40 kHz, which is easily achievable. 

We now turn to the scattering rate T sc associated 
with the fact that dipole-dipole interactions make the 
dark state imperfect on the 421 nm transition. This 
scattering rate is given by T sc = T(Vdd/^dr) 2 , where 
r/27r = 32 MHz is the linewidth of the 421 nm transition 
[49]. To get r sc /27r = 100 mHz and fl 2 r /A = (27r)30 
Hz, we need A/27T = 100 MHz and £ 2 dr/ 27 r = 50 kHz, 
which is easily achievable. Assuming we are tuning to 
the blue of the F' = 17/2 hyperfine level of the excited 
state [see Fig. 1(c)], the resulting detuning from both the 
F' = 15/2 state and the F' = 19/2 state is S ~ (27r)l 
GHz, resulting in scattering from these levels at rates 
~ r(lldr/£) 2 < (27r)100 mHz. 

An important advantage of our proposal is that we are 
not relying on collisions between the atoms, which allows 
us to avoid dipolar relaxation [50]. On the other hand, 
the main limitation of our proposal is the small energy 
scale £fqh/27t = 100 mHz and the resulting stringent 
requirements on the linewidths of the lasers and coher¬ 
ence times of the atoms. An additional requirement is 
to make sure stray magnetic fields are weak enough that 
the Zeeman shifts are small on the scale of dipole-dipole 
interaction strength Vdd/27r = 5 Hz to ensure that the 
crucial jf transitions in Eq. (1) stay resonant. In 
particular, this requires magnetic fields smaller than a 
micro Gauss, which is challenging. At the same time, 
our goal in the present manuscript is to demonstrate the 
controllability and the potential of magnetic dipoles in 
general, and magnetic atoms in particular, for obtaining 
fractional quantum Hall states in general and those in 
C — 2 bands, in particular. We leave it to future work 
to optimize and modify the presented schemes with the 
goal of increasing the interaction energies by working at 
smaller lattice spacing [51-54] and by choosing more op- 
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timal level configurations. Indeed, the maximum achiev¬ 
able interaction for two dysprosium atoms Ai at /2 apart 
is ~ // 0 (10MB) 2 /[47r(Ai a t/2) 3 ] = (27r)70 Hz and occurs 
when they are both in a \J Z \ =8 state. Therefore, opti¬ 
mization of our scheme can likely significantly raise Vdd 
from its current value of (27r)3 Hz. 

To prepare the fractional Chern insulator ground state, 
we envision first tuning ©o and <f>o to the part of the 
phase diagram where the ground state is a superfluid. If 
the superfluid phase exhibits weak quantum correlations, 
then it could, in principle, be approximated by a disen¬ 
tangled state, in which each atom is in a well-defined spin 
state. This can be prepared by adiabatically turning on 
the fields responsible for the tensor shift and the dress¬ 
ing on a time-scale faster than the inverse of Vdd- Next, 
we imagine adiabatically ramping from our prepared low- 
energy-density superfluid to the fractional Chern insula¬ 
tor. This can only be done if the phase transition be¬ 
tween the two states is continuous [55], a question which 
is unknown for the (2,2,1) state. 

The detection of the fractional Chern insulator can be 
carried out by effectively realizing Bragg spectroscopy. 
Specifically, the momentum- and energy-resolved spectral 
function can be measured by attempting to drive the \\) — 
|t) transition with a different spatially dependent Rabi 
frequencies and different detunings. One can then use 
the spectral function to identify the bulk gap and the 
gapless chiral edges [56, 57]. 

VI. OUTLOOK 

Our proposal generalizes naturally to other magnetic 
atoms, such as erbium [58] and, to a lesser degree due 
to a smaller dipole moment, chromium [59]. It can also 
be extended to engineer bilayer fractional Chern insula¬ 
tors in electric dipoles, such as Rydberg atoms and polar 
molecules. 

This proposal also opens up avenues for engineering 
other exotic topological states. In particular, by replac¬ 
ing hopping terms (of the form a\aj + h.c.) with pairing 
terms (of the form at at + h.c.) along the rungs of a fi¬ 
nite ladder, one expects to get parafermionic zero modes 
at the ends of the ladder, by analogy with proposals in 
solid-state heterostructures [60-63] and cold atoms [64]. 
The local replacement of pairing with hopping can be 


engineered by dressing the dipoles along the legs of the 
finite ladder in a way that is different from the rest of 
the system. An introduction of an array of such ladders 
into the (2,2,1) state may then lead, by analogy with 
Ref. [65], to a gapped ground state supporting Fibonacci 
anyons, which are universal for topological quantum com¬ 
puting via braiding alone [1]. In fact, by extending the 
conjecture of Refs. [65, 66], partially supported by nu¬ 
merics [67, 68] (although see Refs. [69, 70]), one might 
expect that even the introduction of uniform pairing on 
top of the (2,2,1) state may give rise to the Fibonacci 
phase. Instead of introducing pairing, one might also be 
able to obtain the Fibonacci state by introducing uniform 
tunneling between two fractional Chern insulator layers 
[66, 68], something we can also engineer in the dipolar 
approach. The use of uniform pairing or tunneling is 
an exciting prospect as the corresponding dressing would 
be significantly simpler relative to the approach involv¬ 
ing an array of ladders [65] and would therefore consti¬ 
tute a particularly promising approach towards universal 
topological quantum computing in ultracold atomic and 
molecular systems. 
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Appendix A: Details behind the C — 1 topological 
fiat band and the v — 1/2 Laughlin state 

In this Appendix, we present the details behind the 
C = 1 topological flat band and the v = 1/2 Laughlin 
state discussed in Sec. II. 

Assuming the tensor Stark shift V t is greater than the 
typical strength Vdd of dipole-dipole interactions, we can 
project the Hamiltonian in Eq. (1) on the 36-dimensional 
Hilbert space, where each atom is in one of the states |1) 
through 16) defined in Fig. 1(c), to obtain 


Hij = (1 - 3cos 2 0) L 2 (n 4 - h 3 )(h 4 - n 3 ) + /&(|15> (51| + |26) (62| + |16) (52| + |25) <61| + h.c.) 

-i/4(|13> <311 + 124) (421 + h.c.) - ±/4(|53> (35| + |64) (46| + h.c.) 
sin 2 9 [e- 2i q M 2 3 (|32> (14| + |23> (41|) - /4(|36) <54| + |63) <45|)} + h.c.} , (Al) 

I 

where /i 4 = (4| J z \4) = |§§, ^26 = (2| J Z \Q) = > Mi 3 = <1|^ + |3) = and /x 53 = (5|J+|3) = 



Assuming the separation r /A between the dark 
states |t), |t) and the bright states is larger than the 
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interaction strength Vdd, we can further project the re¬ 
sulting Hamiltonian on the 4-dimensional Hilbert space 
spanned by |t) = | 0 ) and |t) = o) | 0 ) on each atom: 


Hij = (1 — 3 cos 2 



0)[nihj {/4(k| 2 - \si\ 2 ){\vj\ 2 - \ Sj \ 2 ) - Ifilaiw-Wjisis* + v iV *) + h.c.)} 

+ K w j - l^i 3 ( s * s i + v * v j)) + h c -} 

e~ 2l<t> 11 iil 3 hihj(s*w*VjWi + s*w*ViWj) - 3 {a\ajVjS* + a^ViS*)} + h.c. 


(A 2 ) 


where fij = atm,. 


In Fig. 2, we use @ 0 = 2.35 and <F 0 = 4.13 as 
the angles determining the direction of the DC elec¬ 
tric field. The remaining parameters on the two 
sublattices (which we denote by 1 and 2 in the sub¬ 
script) are parametrized as s 1 / 2 = sin(aq/ 2 ) sin( 6 > 1 / 2 ), 

vi/2 = sin(o: 1/ 2)cos(6» 1/ 2)e^ 1 /2 ) Wl/2 = 

cos(ai/ 2 )e * 71 / 2 with { 6 »i, 0 2 , </>i, </> 2 , «i, 02 , 71 , 72 } = 

{2.28,1.59, -0.05,1.51,0.46,0.11,1.60,1.23}. As ex- 
plained in the main text, we make a further modification 
by changing the sign of w\ and w 2 every other row. 
While this modification does not increase the size of the 
unit cell, it is an important ingredient allowing us to 
obtain a sufficiently flat topological band. 


In passing, we note that we can get the same band- 
structure as for polar molecules [13] in a small electric 
field by adding two additional levels to the linear super¬ 
position composing |t). 


Appendix B: Details behind the C — 2 topological 
fiat band, the Halperin (2,2,1) state, and the 
corresponding phase diagram 

In this Appendix, we present the details behind the 
C = 2 topological flat band and the Halperin ( 2 , 2 , 1 ) 
state discussed in Sec. IV. 

Instead of using the simplest possible level structure for 
obtaining a flat C = 2 topological band with 161 Dy, we 
instead choose to demonstrate the full power and the full 
tunability of dipolar interactions between dressed states 
by realizing a Hamiltonian describing two-component 
hardcore bosons, each of which obeys the same Hamil¬ 
tonian as hardcore bosons in C = 1 bands realized with 
polar molecules [13]. This derivation allows us to estab¬ 
lish a certain degree of equivalence between implemen¬ 
tations of topological flat bands and fractional quantum 
Hall states with different dipoles, both electric and mag¬ 
netic. We will therefore consider the level structure and 
the dressed states shown in Fig. Al, where we assume p 
and q are real, positive, and site-independent. 

Transition matrix elements between |ff) and |t) vanish 
because of the F —> F =b 0,1 selection rule ensuring that 
absence of terms exchanging the two types of excitations 
| ft) and |t). Starting with Eq. ( 1 ), and keeping only those 
terms that will contribute to the interaction involving 
states |t), |t) 5 and |ff), we obtain 


Hij — (1 — 3 cos 2 0) |^(/^6(^6 — ^ 5 ) + M7(f^7 — n%) + M 14(^14 — 7113) + /il6(^16 — ^ 15 )) X 
X (p6(fle — 715) + 127(717 — fig) + ( 1 14(7114 — ff'is) + M 15(^15 — 7116)) 

+/4(|19> (911 + 12,10) (10,2| + |1,10) (921 + |29) (10,1| + h.c.) 

+A*?,n(|l> 11) (11,1| + |2,12) (12, 2| + |1,12) (11,2| + |2,11) (12,1| + h.c.) 

-3/4(135) (531 + |46) (64| + h.c.) - |/4,(|38) (83| + |47) (74| + h.c.) 

sin 2 9 [e- 2i q/4(|36) (54| + |63) (45|) + m| 3 (|84) (37| + |48) (73|)} + h.c] (Al) 

= (1 — 3cos 2 6 ) + }a\a j (iJ,l 9 w*w j p 2 - + v*Vj)q 2 ) + h.c.} 

+ {WhinW*W jP 2 - §/&(StSj + V*Vj)q 2 ) + h.c.}] 


-f sin 2 0 


e 2llp { - \nl 5 {a\ajV*Sj + aia]v*Si)q 2 - lulzfijbjV*Sj + bib]V*Si)q 2 } + h.c. 


(A 2 ) 
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[ 20 } 


iM) _ 

1H} _ 117) 


J5). |3) 


HD 

ID 


I 19 ) F = 15/2 

Hi) F = 17/2 
F = 19/2 

118) _ ill} F = 21/2 


Hi — — — 

_ _ J2L ^ 

— J£L _ 110} —- 

M =-3/2 -1/2 1/2 3/2 5/2 7/2 9/2 11/2 13/2 15/2 17/2 19/2 

u>=iWW 3> + |4> 


>/2 




It) = c |17) J- |18) + a|13) + 6|14) - s|5) + «|6) + J 9) ^ 

ID-) = C |19) + A|15) + F|16) - S\7) + H|8) + 

v 2 v 2 

FIG. Al. The level structure of 161 Dy relevant for implementing the C — 2 flat band and the resulting Halperin (2,2,1) state. 


where the second expression is obtained by making a 
further projection onto \\) = |0), |t) = o) |0), and 
lit) = |0). Here hi = afa;, m; = b\bi, = 

Me(k | 2 - |si| 2 ) + Mi4(|D| 2 - k| 2 ), and fi m>i = p 7 (\Si \ 2 - 
\V i \ 2 ) + p l 5 {\A i \ 2 -\B i \ 2 ). 

To reduce the Hamiltonian governing a and b to two 
copies of hardcore bosons obtained from polar molecules 
in near-zero electric field [13], we define siq = rxi , viq = 

ryi, Wip = T- 2 j|/W( 2 Mi 9 )|, s iQ = IM35/M83I rXi, V^q = 

|/W/h33 \rYi, WiP = rZi\jj, 35 /(2fi ltll )\. We keep |sj| 2 + 
|D| 2 + \u>i\ 2 = 1 and |5j| 2 + |V/| 2 + |Ht| 2 = 1. r is 


a positive real number, c and C are used to keep the 
states normalized. To mimic polar molecules, which are 
electric dipoles and thus have no induced dipole moments 
at small electric fields, we would like to choose a^, 
and Bi such that /i n ^ = /i m ^ = 0. To do this, we reduce 
r from 1 and below until we find p such that there is 
enough population in a$, bi , Ai, and Bi levels to cancel 
dipole moments and for a and b, respectively, 
for both sublattices. For example, for the parameters 
used to produce Fig. 5, r = 0.265 and p = 0.953. 

The resulting Hamiltonian is: 


r 2 /4" tJ 


Hij = (1 — 3 cos 2 6) lajajiwtwj - + v*vj)) + b^W*^ - l(S*Sj + V?Vj )) 


+ | sin 2 0 cos(20) a\dj(v*Sj + s*Vj) + b\bj(V*Sj + S*Vj) 
— | sin 2 Qsm(2(/))i a\dj(v*Sj — s*Vj ) + b\bj (V*Sj — S*Vj ) 


+ h.c.. 


(A3) 


As desired, it has no density-density interactions, and 
both a and b obey the hopping Hamiltonian derived from 
polar molecules [13]. 

In Fig. 5, we use ©o = 0.68 and <Fo = 5.83 as 
the angles determining the direction of the DC elec¬ 
tric field. The remaining parameters on the two sub¬ 
lattices (which we denote by 1 and 2 in the subscript) 
are parametrized as s 1 / 2 = *§ 1/2 = sin(A 1 / 2 ) sin(6 > 1 / 2 ), 


V 1/2 = V 1/2 = sin(a 1/2 )cos(0 1/2 )e^ 1 / 2 , w 1/2 = W 1/2 = 
cos(di/ 2 )e * 7l /2 with { 61 , d 2 , />i, fa, di, a 2 , 71 . 72 } = 

{0.53,0.97,1.36,3.49,2.84,2.03,4.26,3.84}. As in the 
C — 1 example, an additional minus sign is imposed on 
w 1 / 2 and W 1/2 on every other row. It is important to 

repeat that, while a and b here obey the same Hamilto¬ 
nian, one should be shifted relative to the other by one 
lattice site in the X direction. The parameters used in 
Fig. 6 are the same as those in Fig. 5, except ©o and <f>o 
are varied. 
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